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ABSTRACT 



The free streaming of warm dark matter particles dampens the fluctuation spec- 
trum, flattens the mass function of haloes and imprints a fine grained phase density 
limit for dark matter structures. The phase space density limit is expected to imprint 
a constant density core at the halo center on the contrary to what happens for cold 
dark matter. We explore these effects using high resolution simulations of structure 
formation in different warm dark matter scenarios. We find that the size of the core we 
obtain in simulated haloes is in good agreement with theoretical expectations based 
on Liouville's theorem. However, our simulations show that in order to create a signif- 
icant core, (r c ~ f kpc), in a dwarf galaxy (M ~ 1O 1O M0), a thermal candidate with a 
mass as low as O.f keV is required. This would fully prevent the formation of the dwarf 
galaxy in the first place. For candidates satisfying large scale structure constrains (m„ 
larger than ps 1 — 2 keV) the expected size of the core is of the order of 40 (80) pc 
for a dark matter halo with a mass of f0 10 (fO 8 ) M@. We conclude that "standard" 
warm dark matter is not viable solution for explaining the presence of cored density 
profiles in low mass galaxies. 
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1 INTRODUCTION 

The formation of structure in the universe is driven by 
the mysterious dark matter component whose nature is 
still unknown. Over the last decades, the hierarchical cold 
dark matter model (CDM) has become the standard de- 
scription for the formation of cosmic structures. It is 
in excellent agreement with recent observations, such as 
measurements of th e cosmic microwave b ackground and 
large scale surveys (|Tegmark et al.l 120061 ; iKomatsu et al.l 
l201ll ). However, there are a number of inconsistencies on 
sub-galactic scales that arise within the CDM scenario. 
Firstly, the amount of substructure in Milky Way sized 
haloes is overpredicted by roughly one order of magnitude 
l|Klvpin et alj 11999b iMoo're et al.lfl999T ). Secondly, the cen- 
tral de nsities of CDM haloes in simulations show a cuspy be- 
havior (Moore 1994; Florcs & Primack 1994; Dicmand et al. 
120051 ; iMaccio et all 120071 ; ISpringel et ail |200St ). whereas 
the density profiles inferred from ga laxy rotation curves 
point to a core like structure (e.g. Ide Blok et al.l l200ll ; 
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iKuzio de Narav et al.l 12009b lOh et al.l l201ll). Furthe r more , 



recent studies |Tikhonov et al. 20091 ; Zavala et al.l 12009b 
iPeebles fc Nusserll2010f ) re-emphasized that also the popula- 
tion of dwarf galaxies within voids is in strong contradiction 
with CDM predictions. 

One possible solution to these issues is that the dark 
matter particle is a thermal relic with a mass of order one 
keV. The most prominent representatives of such warm dark 
matter ( WDM) candidates are the sterile neutrino and the 
gravit ino (|Abazaiian fc Koushiappasl 120061 ; iBovarskv et al.l 
2009a|), whose presence is also motivated by particle the- 



ory (e.g. Dodelson & Widrow 



Takav ama fc Yamaguchl koOO) 



1994; iBuchmiiller et alj|2007l : 



)W 



Non-zero thermal velocities for WDM particles lead 
to a strong suppression of the linear matter power spec- 
trum o n galactic and s ub -galac tic scales ( Bond et al 
19801 ; iPagels fc Primackl 1 19821; [Dodelson fc Widro 
1994 iHogan fc Dalcantonl l2000b IZentner fc Bullock! [20031 ; 
Abazaiian 20061 ; IViel et al.ll2005l ). and erase all primordial 
density perturbations smaller than their free-streaming 
scale A f s . Below this scale no structure is expected to form, 
at least not in the usual bottom-up scenario. However, the 
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effective suppression of halo formation already happens well 
above Xf a and is entirely describe d by the WDM particle 
mass (see lSmith fe Markovidl201ll . and references therein). 

Recent observational constraints coming from X- 
ray background measurements and Ly-a forest analy- 
sis set the allowed mass interval r oughly between 2 
and 50 keV (e.g. IViel et alJ 120051; ISeliak et al.l l200d 
lAbazaiian fc Koushiappasl 12009 ; iBovarskv et al.l l2009bl lcf) I 1 
As a complementary study Maccio & Fontanot (2010, see 
also Polisensky & Ricotti 2011) compared the subhalo abun- 
dance of a Milky Way like object in different numerical warm 
dark matter realizations with observed satellite galaxies re- 
ported by the SDSS data and set a lower bound for a ther- 
malized particle of itiwdm > 2 keV. 

Another important characteristic of a WDM scenario is 
the possibility to naturally obtain cored matter density pro- 
files. According to Liouville's theorem for collisionless sys- 
tems, the finite phase space density of the cosmic fluid stays 
constant through cosmic history. In WDM the dark mat- 
ter fluid is described by a Fermi-Dirac distribution, whose 
absolute value is fixed at the time of decoupling when the 
fluid becomes collisionless. Structure formation then hap- 
pens through a complex process of distortion and folding of 
the phase space sheet. Since it is not possible to measure this 
fine-grained phase space density in simulations, one usually 
defines a coarse-grained or pseudo phase space density (e.g. 
iTavlor fc Navarroll200ll ) 

(i) 

where p is the mean density and a is the one-dimensional 
velocity dispersion within some small patch of the simu- 
lation. The quantity Q corresponds to an average density 
of a small (but not microscopic) phase space volume and 
is not constant anymore. However, because of the way the 
phase space sheet is distorted, the value of Q can only de- 
crease during structure formati on and will never exceed it s 
initial value set at decoupling l|Dalcanton fc Hoganll200~it ). 
This fundamental upper limit also holds for the local phase 
space density within virialised haloes at redshift zero and 
has a direct consequence on the density profile in real space. 
Since the velocity dispersion does not grow in the inner part 
of a halo, the real space density profile must become con- 
stant with a core size dep ending on the specific WDM model 
jTremaine fc GunrJ[l979T ) . 

Due to this effect of core formation, the WDM sce- 
nario has been suggested as a solution to the long stand- 
ing core-cusp problem of dwarf galaxies. In fact, obser- 
vational measurements favor cored dark matter profiles 
in low surface brigh t ness galaxies within the loca l grou p 
jSalucci et all l201ll ; iKuzio de Narav fc Kaufmannl l201lf) . 
However, previous theoretical/analytical studies argue that 
the cores produced by warm da rk matter might be too small 
to explain the observations (|Kuzio de Narav et al.l l20ld : 
IVillaescusa-Navarro fc Dalalll201lj ). 

Numerical N-body simulations have been used to 



better understand the properties of virial ized objects in 



the Warm Dark Matter scenarios (e.g. Bode et al.l 


2001 


Kncbe et al. 2003; Wang & White 200 


7l; IZavala et al.l 


2009 


iTikhonov et al. 


120091; ISchneider etaL 


l201ll). High resolu- 



tion simulations of single objects have studied the suppres- 
sion of t he galactic satellit e formation due to free stream- 
ing (e.g IColrn et al l |2000|; iGotz fc Sommer-Larsenl |2002| ; 



iKnebe et al.l 120081 ; iMaccio fc FontanotJ l20ld ). in order to 
reconcile the observed dwarf galaxy abundance with the 
prediction f r om D ark Matter based theories. More recently 
IColfn et alJ (|2008h used N-body simulations to study the ef- 
fects of primordial (thermal) velocities on the inner structure 
of dark matter haloes, with particular attention on the for- 
mation of a possible central density core. They used thermal 
velocities of the order of 0.1 and 0.3 km/s, without linking 
them to any particular WDM model, since the aim of their 
work was to explore the general effect of relic velocities of 
the DM structure. Unfortunately their combination of res- 
olution and choice for relic velocities was not sufficient to 
directly test simulation results against core radii predicted 
by phase-space constraints. 

In this work we want to extend and improve on these 
previous studies. We will use high resolution N-body simula- 
tions to explore the sizes of density cores in WDM and their 
dependence on the WDM candidate masfl We will explore 
several models for WDM ranging from 2 to 0.05 keV. We 
will consider separately the effects of a warm dark matter 
candidate on the power spectrum and on the relic velocities, 
trying to disentangle the various consequences of these two 
different components. Our higher numerical resolution will 
allow us to directly see the formation of a density core, with 
a size well above the numerical resolution for the warmer 
candidates. We will then revise the theoretical arguments 
for the formation of cored profiles in WDM and perform a 
direct comparison between the core sizes in our simulations 
and the ones predicted from phase space constraints. 

The paper is structured as follows: In section 2 we dis- 
cuss the setup of our simulations and the way we implement 
thermal velocities. Section 3 is dedicated to the presenta- 
tion of our results in terms of the phase space limit and its 
influence on the density profile of dark matter haloes. A con- 
clusion and a summary of our work is finally given in section 
3. 



2 SIMULATIONS 

Numerical simulations have been carried out using pkd- 
GRAV, a treecode w ritten by Joachim Stadel and Thomas 
Quinn (|Stade]||200ll V The initial conditions are generated 
with the GRAFIC2 package (|Bertsch ingcr 200j]). All simula- 
tions start at redshift Zi — 99 in order to ensure a proper 
treatment of the non linear growth of cosmic structures. 

The cosmological parameters are set as follows: 
fi A =0.727, fi m =0.273, n b =0.0U, h = 0.7 and cr 8 = 0.8, 



1 In some of these analysis the warm dark matter particle is as- 
sumcd to be a resonantly produced sterile neutrino llShi fc Fuller! 
1999). We have converted these mass limits into limits for a fully 
thermali zed particle , such as the gravitino, using the formula pro- 
vided bv lViel et all \200$) . 



2 In the present work we only considered a very simple WDM 
model; it is worth commenting that there are more complex 
and physically motivated models discussed in the literature (e.g. 
warm+cold dark matter. i B oyarsky et al.| l|2009d|) or co mposite 
dark matter iKhlopovl ^OOd^lKhlopov fc Kouvarij teuOSl ')') 
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Figure 1. Density map of the large scale (low resolution) simulations (L=40 Mpc) at redshift zero. From left to right: CDM, and two 
WDM with a cut in the power spectrum for a m v mass of 0.2 and 0.05 keV respectively. These last two simulations have not been used 
in this paper and are presented only for illustration purposes, see section b.ll for more information. 



and are in good agreement w ith the recent WMAP mission 
results l|Komatsu et al.|[201lh . 

We start by running large scale simulations of a cos- 
mological cube of side 40 Mpc, using 2 x 256 3 dark matter 
particles. This was done for two different models, a stan- 
dard LCDM and a Warm Dark Matter model with a warm 
candidate of mass 2 keV produced in thermal equilibrium. 

To compute the transfer functio n for WDM models we 
used the fitting formula suggested by iBode et alJ (|200lh : 



Table 1. Simulations parameters 



Label 


m v 




iV vir 






(keV) 


(keV) 


(106) 


(10 12 M Q ) 


CDM 


oo 




10.2 


1.42 


WDM1 


2.0 


2.0 


8.6 


1.22 


WDM2 


2.0 


0.5 


8.4 


1.20 


WDM3 


2.0 


0.2 


8.5 


1.21 


WDM4 


2.0 


0.1 


6.7 


0.93 


WDM5 


2.0 


0.05 


4.9 


0.71 



pWDM 

T2 ( k ) = J^M =[l + (akfT W/ " (2) 

where a, the scale of the break, is a fu nction of the WDM pa- 
ram eters, while the inde x v is fixed. IViel et all (|2005l ) (see 
also lHansen et al.l l2002), using a Boltzmann code simula- 
tion, found that v — 1.12 is the best fit for k < 5 h Mpc - , 
and they obtained the following expression for a: 

(3) 

We used the expression given in Eq. ([3]) for the damping of 
the power-spectrum for simplicity and generality. More ac- 
cura te expressions for the damping of sterile neutrinos exist 
(e.g. lAbazaiianll2006h and show that the damping depends 
on the detailed physics of the early universe in a rather 
non-trivial way. The initial conditions for two simulations 
have been created using the same random phases, this means 
that there is a one-to-one correspondence between CDM and 
WDM haloes. 

We then select one candidate halo with a mass similar 
to our Galaxy (M ~ 1O 12 M0) and re-simulated it at higher 
resolution. These high resolution runs are 8 3 times more re- 
solved in mass than the initial ones: the dark matter particle 
mass is m p — 1.38 x 10 Mq, where each dark matter particle 
has a spline gravitation softening of 355 pc. This single halo 
has been re-simulated in several different models, all simula- 
tions are summarized in Table[T]and three of the simulations 
are shown in Figure 1. 



2.1 Streaming velocities 

Particles that decouple whilst being relativistic are expected 
to retain a thermal velocity component. This velocity can be 
computed as a function of the WDM candidate mass ( m v ) 
according to the following expression (|Bode et al.ll200lh : 

v (z) m ~f^y f h y (1.5^ fkeV\% _ a 

— = - 012 U^ U^J UJ UrJ kms 

(4) 

where z is the redshift. The distribution functio n is (e p ^ T " + 
l)" 1 until the gravitational clustering begins (|Bode et al.l 
l200ll ). 

This formalism is correct for the "real" dark matter el- 
ementary particles (e.g. a sterile neutrino). In the N-body 
approach we use macro particles (with masses of the order 
of 10 j Mq) to describe the density field. These macro parti- 
cles effectively model a very large number of micro particles. 
Given that the velocities described in Eq. Q have a random 
direction the total velocity of the macro (Nbody) particles 
should effectively be zero. Hence, it is not correct to directly 
use Eq. @ to assign "thermal" velocities to simulation par- 
ticles. 

On the other hand, the net effect of the thermal ve- 
locities is to create a finite upper limit in the phase-space 
distribution (PSD) due to their initial velocity dispersion 
(a). What we are interested in is to recreate the same PSD 
limit in our simulation, and then study its effects on the 
dark matter halo density distribution. In order to achieve 
this goal we proceed in the following way. From Eq. ([4]) we 
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compute the rms velocity: a(z) = 3.571«o(z), we then create 
a Gaussian distribution centered on zero and with the same 
rms a. Finally we randomly generate particle velocities from 
this distribution and assign them to our macro particles. 

As detailed in section 13.11 there is a direct connection 
between m„ and the expected size of the dark matter dis- 
tribution core. This core is only due to the presence of ther- 
mal velocities and not, in first approximation, to the cut 
in the power spectrum described by Eq. (2j. Cutting the 
power spectrum changes the merger history of the dark mat- 
ter halo but does n ot affect the density profile significantly 
l|Moore et al.ll 19991 ). This implies that in order to study the 
effect of different values of m v (and hence Vo) it is sufficient 
to "play" with Eq. Q leaving all other simulation param- 
eters unaltered. Following this approach we have generated 
several simulations using the same cut in the power spec- 
trum {m v ) but different initial thermal velocities (m„ jV ci), 
as detailed in Table 1. 



3 RESULTS 

Density profiles for the cold dark matter run and the four 
warm dark matter realizations (WDM1-WDM5) are shown 
in Fig. [2] The profiles show a monotonic decrease of the cen- 
tral density as a function of the the temperature of the dark 
matter can didate. Cold candidate show the usual cuspy be- 
havior (e.g. iDubinski fc Carlberslll99lf) . while warmer can- 
didates present a lower central density that becomes a clear 
core for m„ jVC i = 0.05 keV, with a size of several kpc. 

Fig. [3] shows the time evolution of the density profile in 
the WDM5 simulation. The core is already in place at high 
redshift z = 1.6, and does not evolve substantially until 
z = 0. The profile only changes at large radii (r > 50) kpc, 
as a consequence of the assembly of the external part of the 



halo that is consistent with a typical CDM halo in the outer 
regions. 

As already mentioned the theoretical explanation for 
the formation of a core is related to the presence of a max- 
imum in the phase space density distribution. This maxi- 
mum is clearly visible in Fig. [4] where we plot the phase 
space density proxy Q = p/a 3 , for three different models, 
namely CDM, WDM3 and WDM5. For this latter model 
the phase-space proxy Q shows a large core that extends 
about 10 kpc. The WDM3 model also shows a strong flat- 
tening of the Q profile, consistent with a core distribution. 
On the other hand the CDM phase-space distribution is well 
fitted by a single power law pr ofile on the whole ra nge, in 
agreement with previous results ISchmidt et all |2008l ). 

Figure [S] shows the time evolution of the PSD for our 
warmest candidate (i.e. thermal velocities for a 0.05 keV 
mass particle). The solid (blue) line shows the Q radial pro- 
file in the initial conditions [z = 99). This value has been 
calculated using only high resolution particles that end up 
within 1.5 times the virial radius of the halo at z = 0. These 
particles define a lagrangian region with a comoving volume 
of about (6Mpc) 3 . The other (red) lines represent the PSD 
profile at different redshifts (from 1.6 to 0) and have been 
computed using all particles within the virial radius of the 
halo. All quantities in the plot are in physical units. The 
phase space distribution shows very weak evolution with al- 
most no evolution at all from z = 99 to z = 1.6. In the same 
plot we also show the theoretical maximum phase-space den- 
sity achievable by this model (see Eq. Q for a rigorous def- 
inition of Qmax)- As we will see later (in section [3TTJ the 
value of this maximum phase space density depends on the 
local value of Q m . The two (black) lines show predictions 
for Qmax for the average matter density in our simulation 
(fl m = 0.273) and for its local value, which we measured by 
dividing the mass contained in the halo lagrangian region 
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Figure 4. Phase space density profile for the CDM, WDM3 and 
WDM5 models at z = 0. 

(defined above) by its volume. This local matter density pa- 
rameter turned out to be Q m> i OC ai ~ 0.0344 0. 

This second theoretical prediciton, that took into ac- 
count the local value of the matter density parameter, is 
in very good agreement with the simulation results. It also 
shows that global prediction on the size of cores in WDM 
could be affected by the local overdensity. We will come back 
to this issue later. 

In order to quantify the flatness (and the core size) of 
WDM profiles we have fitted all our density profiles with 
the following param etric description, originally presented in 
IStadel et all (|2009h : 

p(r) = p e X p(-A[ln(l + r/7? A )] 2 ). (5) 

In this parameterization the density profile is linear 
down to a scale R\ beyond which it approaches the cen- 
tral maximum density po as r — > 0. We also note that if one 
makes a plot of dhip/dln(l + t/Ra) versus ln(l + r/RA) then 
this profile forms an exact straight line with slope 2A. 

This fitting function is extremely flexible and makes 
possible to reproduce at the same time both cuspy profiles 
like the one predicted by the CDM theory, as well as, highly 
cored profiles, like in the WDM5 case (as shown in Fig. [5]). 
The values of the parameter are obtained via a \ 2 minimiza- 
tion procedure using the Levenberg & Marquart method. 
From now on we will use the value of the fitting parameter 
R\ as the fiducial value of the central density core in sim- 
ulated profiles (r cor0iS , hereafter). The r cor0jS values for all 
our haloes are reported in the first column of Table 2. 
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Figure 5. Time evolution of the PSD radial profile for the WDM5 
model. The two (black) lines represent the theoretical prediction 
for the maximum value og Q according to equation [7] The upper 
one is for the global value of the matter density paramter (f2 m ) 
while the lower one is for its local value (Q m i oca i = 0.0344, See 
text for more details). 
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Figure 6. Density profiles for CDM and WDM5 and their fit 
using Eq. 10. 



3 This low value is not surprising since the mass of the halo with 
1.5 R v i r is of the order of 5 X 1O 12 M0 and it is distribute on a 
volume of (6Mpc) 3 in the IC. This is partially a consequence of 
having chosen in our study a very isolated halo 
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3.1 Comparison with theoretical predictions 

In iTremaine fe Gunnl (| 19791 . TG79 hereafter) limits on 
the mass of a neutrino are derived from the maximum 
phase space density of a homogeneous neutrino background, 
with the further assumptions that neutrinos form bound 
structures and that their central regions can be well- 
approximated by an isothermal sphere. 

Assuming a Maxwellian velocity distribution they ob- 
tained the maximum phase space density: 



P 4 
— r oc m v 



(6) 



where m v is the mass of the (warm) dark matter candidate. 
This limit has been then used in several follow up papers 
to estimate the size of density cores in warm dark matt er 
haloes (e.g. iDalcanton fc Hoganll200ll ; IStrigari et alll2006j ) . 

While the calculation from TG79 is formally correct, it 
cannot be directly applied to our numerical simulations (and 
more in general to any numerical simulation). The reason 
is that when deriving Q max , TG79 tried to constrain at the 
same time both the phase space density and the total matter 
density (Q m ) in warm (or hot) dark matter. In doing that 
they correctly assumed Qwdm oc m v . This is not what is 
usually done in cosmological simulations, where the value of 
tt m is kept constant among different models (e.g CDM vs 
WDM) and fixed to match cosmic micro wave based results. 

It is then worthwhile calculating the value of Qmax un- 
der the assumption of a fixed value for the total dark matter 
density. We can start from the definition of Q: 



Vim Per 



(7) 



where p cr = 2.775 x lO 11 /i 2 Af Mpc -3 is the critical density 
of the Universe. Following TG79 we can rewrite the denom- 
inator according to: 



kT v 
cm v 



(8) 



The temperature of the WDM particle (T„) can be computed 
from the temperature of the p hotons if the WD M candidate 
has a mass below 1 MeV (e.g. IWeinberelll972l ). this gives: 



1/3 



T-v = 1.9K 



(9) 



If we now combine Eqs. l[7[l. Q and ((9J we obtain a new 
expression for maximum phase space density: 



Qmax = 2.74 x 10" 



M pc" 



\keVJ VO-25 y V0.7 7 (kms- 1 )" 3 ' 

(10) 

Aside from numerical factors, the most important difference 
with respect to the expression derived by TG79, is the dif- 
ferent scaling of Q with the WDM particle mass, a direct 
consequence of assuming a fixed (WDM mass independent) 
matter density parameter TO . 

Finally the maximum pha se space density can be con - 
verted in a 'core' size following Irlogan fc Dalcantonl (|200d ) : 



2 

^corc,t 



4:nGQmax < al , > 1/2 ' 



(11) 



Table 2. Size of density cores using different methods. See text 
for a more detailed explanation 



Label 


^*core,s 


^*core,Q 


f*core,t 




(kpc) 


(kpc) 


(kpc) 


CDM 


< 0.4 


< 0.4 


CO 


WDM1 


< 0.4 


< 0.4 


0.009 


WDM2 


< 0.4 


< 0.4 


0.068 


WDM3 


0.42 


< 1.1 


0.272 


WDM4 


1.63 


1.80 


0.69 


WDM5 


4.56 


4.85 


2.79 



where Ohaio is the velocity dispersion (i.e. the mass) of the 



simulated dark matter halo. Values of r coro , t for our simu- 
lated haloes are reported in the last column of Table [2] 

In the following we will compare this theoretical value 
of the core (r corc>t ) with two different core sizes than can 
be estimated directly from the simulations. The first one is 
given by the R\ parameter obtained by fitting the numerical 
density profile, as shown in Fig. [5] and we will refer to this 
value £LS Tcore.s- The second one is obtained by computing 
Qmax from the simulated density profile (as shown in Fig. 
2} and then inserting this value in Eq. (|lip . we named this 
second parameter r core ,Q. 

Results for the three definitions of the core size for all 
our simulations are summarized in Table [2] Overall the three 
different estimators for the core size are in fairly good agree- 
ment. r C: Qmax gives on average a larger value for the core, 
for the WDM4 and WDM5 runs, while for the WDM3 sim- 
ulation is only able to give an upper value, since there is not 
a clear indication of convergence towards a maximum value 
in the Qmax profile, as shown in Fig. [4] 

Fig. [7] shows the comparison of the cored found directly 
in simulations (r corCiS ) with the core predicted by the above 
simple theoretical argument (r CO re,t). In the plot there are 
three independent theoretical predictions for the core value. 
The solid and dashed lines are obtained from Eqs. (|ll|l and 
(|10p . using two different values for the matter density pa- 
rameter Q m , namely the global one (0.273) and the local 
one (0.034). Overall numerical results for WDM3, WDM4 
and WDM5 are in very good agreement with the theoreti- 
cal expectations from Eqs. (JTTJ) and (fT0)l . The WDM1 and 
WDM2 simulations only put upper limits on the size of the 
core, since the values of R\ we obtain from fitting the den- 
sity profile fall below the simulation softening (the dashed 
black line in the figure). The plot also shows that using the 
local value of the density parameter gives a better fit to the 
data. On the other hand the core size only depends on the 
square root of Q m , so the difference between the two core 
size predictions is not very large. A third theoretical value 
(dotted line) is determined using the original argument by 
TG79, where the matter density parameter is linked to the 
mass of the warm candidate. In this case the theoretical rela- 
tion seems to be "tilde" with respect the simulation results, 
and it tends to over estimate the size of the density core in 
simulations with a fixed fl m parameter. 

Using our new determination of the core size as a func- 
tion of the warm dark matter mass we compute the expected 
value of r core for t he typical halo mass (5 x 10 8 Mq, see 
iMaccio et al.l ([2010)) of dwarf galaxies orbiting the Milky- 
Way. Results are shown in Fig. [8] the grey shaded area takes 
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m v (keV) 

Figure 7. Comparison between core size in simulations (open 
symbols) and the theoretical expectation for a M = 10 12 AfQ 
halo (solid line). The dashed line is the gravitational softening of 
our simulations. All points below this line should be considered 
as upper limits on the core size. 




0.01 0.1 1 10 



m„ (keV) 

Figure 8. Expected core size for the typical dark matter mass 
of Milky Way satellites as a function of the WDM The 
shaded area takes into account possible different values of the 
local density parameter 0.05 < Q m < 1.4. The vertical dashed 
line shows the current limits on the WDM mass from large scale 
structure observations. 



into account possible different values of the local matter den- 
sity parameter in the range 0.05 < f2 m < 1.4. 

From the figure it is clear that a core of as 1 kpc would 
require a wdm mass of the order of 0.1 keV, well below 
current observational limits from large scales. 

If we assume a warm dark matter particle mass of 
m„ ~ 2 keV (reprensented by the dashed vertical line), 
in agreement w ith several astrophysical constraints (e.g. 
IViel et all 120081 ). the maximum core size we can expect 
ranges from 10 pc for a massive, MW-like halo (see also 
figure [7|, to 30-50 pc for a dwarf galaxy like halo. Finally, 
in predicting the core size for satellite galaxies in the MW 
halo, it must be taken into account that due to stripping and 
tidal forces satellites c an lose significant mass after accreting 
into larger haloes (e.g. lPenarrubia et al . 2008; M accio et al.l 
|2010| ). This implies that the halo mass we may infer today 
for those galaxies is only a lower limit on the mass they had 
before accretion, which is the one to be used (as Ohaio) in 
Eq.ECD 



4 CONCLUSIONS 

We have used high resolution N-body simulations to exam- 
ine the effects of free streaming velocities on halo internal 
structure in warm dark matter models. We find: 

• The finite initial fine grained PSD is a also a maximum 
of the coarse grained PSD, resulting in PSD profiles of WDM 
haloes that are similar to CDM haloes in the outer regions, 
however they flatten towards a constant value in the inner 
regions. This is in agreement with previous studies based 



on simulations (Colm et al.ll2008l) a nd theoretical arguments 
IVillaescusa-Navarro fc Dalall (|201ll ). 

• The finite PSD limit results in a constant density core 
with characteristic size that is in agreement with theoretical 
expectations i.e. following Tremaine & Gunn 1979, especially 
if value of the local matter density is taken into account. 

• The core size we expect for thermal candidates allowed 
by independent constraints on large scales (Lyman-a and 
lensing, m„ w 1 — 2 keV), is of the order of 20-60 pc. This is 
not sufficient to explain the ob served cores in dwarf galax- 
ies that are around kpc scale l|Walker fe Penarrubia|[201ll ; 
IJardel fc Gebharj2012h . 

• Our results show that a core around kpc scale in 
dwarf galaxies, would require a thermal candidate with a 
mass below 0.1 keV, ruled out b y all large scale structure 
constraints JSeliak et al.l l200fj ; Miranda fc Macciol 120071 ; 
IViel et al.ll2008l ). Moreover with such a warm candidate, the 
exponential cut-off of the Power Spectrum would make im- 
possible to obtain these d warf galaxies in the first place (e.g. 
IMaccio fc Fontanotir2010l ). 

• All together these results lead to a nice "Catch 22" 
problem for warm dark matter: If you want a large core you 
won't get the galaxy, if you get the galaxy it won't have a 
large core. 

We conclude that the solution of the cusp/core prob- 
lem in local group galaxies cannot completely reside in sim- 
ple models (thermal candidates) of warm dark matter. If 
cores are required then it seems that baryo nic feedback (e.g. 
iGovernato et al.ll201ol : IMaccio et alll2012h is still the most 
likely way to alter the density profile of dark matter and 
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hence reconcile observations with cold/warm dark matter 
predictions. 
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